EPJ manuscript No. 

(will be inserted by the editor) 



OO 
O 

o 

(N 
> 

o 
in 



Can the RHIC J/ip puzzle(s) be settled at LHC? 

L. Bravina 1 , A. Capella 2 , E.G. Ferreiro 3 , A.B. Kaidalov 4 , K. Tywoniuk 1 , and E. Zabrodin 1 

Department of Physics, University of Oslo 
0318 Oslo, Norway 

Laboratoire de Physique Theorique a , Universite de Paris XI, Batiment 210, 
91405 Orsay Cedex, France 

Departamento de Fisica de Partfculas, Universidad de Santiago de Compostela, 
15782 Santiago de Compostela, Spain 
Institute of Theoretical and Experimental Physics 
RU- 11 7259 Moscow, Russia 



43 

Oh 



> 
O 
On 

o 



OO 

o 



X 



Received: date / Revised version: date 

Abstract. One observes strong suppression effects for hard probes, e.g. the production of J/tj) or high-pT 
particles, in nucleus-nucleus (AA) collisions at RHIC. Surprisingly, the magnitude of the suppression is 
quite similar to that at SPS. In order to establish whether these features arise due to the presence of a 
thermalized system of quarks and gluons formed in the course of the collision, one should investigate the 
impact of suppression mechanisms which do not explicitly involve such a state. We calculate shadowing 
for gluons in the Glauber-Gribov theory and propose a model invoking a rapidity-dependent absorptive 
mechanism motivated by energy-momentum conservation effects. Furthermore, final state suppression due 
to interaction with co-moving matter (hadronic or pre-hadronic) has been shown to describe data at SPS. 
We extend this model by including the backward reaction channel, i.e. recombination of open charm, 
which is estimated directly from pp data at RHIC. Strong suppression of charmonium both in pA and 
AA collisions at LHC is predicted. This is in stark contrast with the predictions of models assuming QGP 
formation and thermalization of heavy quarks. 

PACS. 13.85.-t Hadron-induced high- and super-high-energy interactions (energy > 10 GeV) - 25.75.- 
q Relativistic heavy-ion collisions - 25.75.Cj Photon, lepton, and heavy quark production in heavy ion 
collisions 



5-j | 1 Introduction The RHIC results on the centrality-dependent nuclear 

modification factor of charmonium in Au+Au collisions at 

y/s — 200 GeV T contained two puzzling features. The 

Charmonium production off nuclei is one of the most pro- firgt wag the fac( . that the suppression at m i d -rapidity co- 

mismg probes for studying properties of matter created m indded ^ the ^ of suppression in Pb+Pb co i Mo ns 

ultrarelativistic heavy ion collisions. Being a heavy par- ftt ^ = u 3 GeV for the same number f participants 0. 

tide, it can be used as a probe of the properties of the gince the medium produced at RHIC is denser and lives 

medium created in these collisions, such as the intensity longer than &t gps Qne WQuld; a priori; expect & stronger 

of interactions and possible thermalization. suppression with increasing energy. The second puzzling 

The RHIC era m charmonium physics has brought f eature was the stronger suppression at forward rapidity 

to light further aspects beyond what was known during than at y = 0. Once again, the medium in the central part 

the SPS experiments. The gold plated signal of charmo- of the collision is the most dense and from simple argu _ 

nium suppression due to colour screening in a quark-gluon ments we would expect an opposite behaviour than what 

plasma was shown to be a characteristic of models which is seen in the data Additionally, data on charmonium in 

did not assume such a state and recently the study of d+Au co n isions at the same energy revealed a decrease of 

charmonium regeneration, or recombination, has become the suppreS sion due to cold nuclear matter compared to 

a field of active study. In this sense the results from RHIC lower energies 
establish an important step, bridging the vast energy jump 

to the future LHC experiments. These seemingly puzzling features in the data bring to 

light two new aspects of charmonium physics in nuclear 

a Unite Mixte de Recherche UMR n° 8627 - CNRS reactions, namely the role of initial state effects and the 

Correspondence to: konrad.tywoniuk@fys.uio.no role of possible secondary J/ljj production from recombi- 
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nation of open charm. At RHIC these effects are visible 
and important for the detailed description, at LHC their 
magnitude, especially the gluon shadowing, will be much 
larger. We argue that only the full description of charmo- 
nium dynamics from SPS to LHC energies will shed new 
light on the underlying physics of dense partonic systems. 



2 Baseline: initial state effects 

In order to quantify the modifications of charmonium due 
to the presence of a dense partonic medium, one usually 
starts from the simpler proton- nucleus (pA) collisions. At 
least two modifications in the pA case compared to pp are 
known to come into play. On the one hand, charmonium 
suppression is seen to scale with L, the traversed path 
length by the initially formed cc pair within the nuclear 
medium. In the Glauber model, this leads to a suppres- 
sion factor oc exp (— pLa a bs), controlled by the so-called 
absorptive cross section, a a bs ■ In P+Pb collisions at ^fs = 
19 GeV this cross section was found to be ~ 4.5 mb [3]. On 
the other hand, the nuclear gluon distribution is modified 
at high energies, leading to 
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at values of Bjorken-x < 0.1, called gluon shadowing. Usu- 
ally one assumes a gluon fusion mechanism for the char- 
monium production, and therefore this effect is of crucial 
importance. 

Recently, it has been realized that the nuclear suppres- 
sion of charmonium in pA collision exhibits a non-trivial 
energy dependence (see e.g. [HE]) namely, that <j a b s is 
decreasing with energy , 17.3 < y/s < 45 GeV, in con- 
trast to elementary theoretical expectations [5J. At RHIC, 
the situation is more complex due to the presence of ad- 
ditional gluon shadowing. The absorptive cross section 
introduced above is replaced by a rapidity-independent 
(^break-up which takes values in the range {0.7,4.5} mb, 
depending on which shadowing parameterization is chosen 
[318] . This approach seems to be in conflict with the strong 
^-dependence of J/tjj suppression observed in [9]J10| . al- 
though the accessible range of xp at RHIC is quite limited. 
In summary, there are still many open questions which 
seems to lack an explaination within a unified framework. 

The observed trends are, in fact, in line with expec- 
tations from the Glauber-Gribov theory [Tl"l|12| . In this 
approach the quantity a a b s controls the contribution of a 
certain set of diagrams valid at low energies, which corre- 
sponds to the Glauber model. More involved diagrams be- 
come, of course, important at high energies [13j . They cor- 
respond to a space-time picture where the initially small 
projectile evolves into a large fluctuation which can inter- 
act coherently with all of the constituents of the target 
nucleus. This leads to nuclear shadowing. The transition 
from the former, planar, to the latter, non-planar, regime 
is governed by a critical energy scale 
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where M c5 is the mass and x+ = (y^xp + AM%/ s + xp)/2 
is the longitudinal momentum fraction of the heavy sys- 
tem. This transition signals the breakdown of the semiclas- 
sical probabilistic picture of longitudinally ordered multi- 
ple scattering as one goes to high energies, s > sm- In 
this scenario, the diagrams that are controlled by a a bs are 
cancelled in the summation and, instead, shadowing of 
nuclear partons appear [14] . 

At finite energies there are also well know corrections 
due to conservation of energy- momentum 15, 11 . These 
are, of course, most prominent when xp — > 1, but can also 
lead to suppression at mid-rapidity. 

These features lead us to propose the following model 
for charmonium suppression at all energies, see |16j for the 
details. The suppression of charmonium in pA compared 
to A elementary pp collisions is given by 

Rf A {xF,s) = i / d 2 bS FE (x + ,b)S shad (x 2 ,Q 2 ,b) . (3) 



S FE (x+,b), is 



The finite energy (FE) suppression factor, 
given by [12] 

<J fe (t h\ - I I t 1 ~ exp(-&TA(&))] for s < s M 
1*+,*;-^ Ta ^ exp (_gr A (&)/2) for s > 



> s M 
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where the nuclear profile function, Ta(6), is normalized 
to the atomic number, A, and the shadowing correction, 
S shad (x2,Q 2 ,b), for the gluons in the nucleus carrying a 
momentum fraction xi is taken from [17] . Note, that no 
antishadowing effects are included in this model. 

The "generalized" absorptive cross section that enters 
eq. ([4} is a function of x+ [TT1[16] 



a(x+) = [(1 - e)<P(t min ) + , 



(5) 



where a c5 is the total cc — N cross section, 7 = 2 for 
charmonium [llj and e is a parameter characterizing the 
amount of longitudinal momentum lost in each rescatter- 
ing. The "form factor" 



${t m in) = exp{-(a^/a;2) 2 } , 



(0) 



where x\ « 0.1 corresponds to the critical energy scale in 
eq. @ , controls the transition from the planar to the non- 
planar regime. Note, that in the low-energy limit and at 
mid-rapidity, where no shadowing corrections exist, S shad = 
1, the ordinary absorption in this model is given by 



(7 ► 



[(1 



e) + eM c 2 5 / S ] a c 
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whereas the Glauber result, without energy- momentum 
conservation, would be simply a a bs — (1 — e)cce- In the 
high-energy limit at mid-rapidity , on the other hand, 



a — ► , 



(8) 



so that S FE — 1 and the suppression is fully governed by 
the shadowing corrections. The parameters cr c5 and e have 
been extracted from low energy data in [16j . 
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Fig. 1. Rapidity dependence of J/ip suppression in pA colli- 
sions at RHIC and LHC. Data are from 08]. 

In figure [T] we compare our calculations of gluon shad- 
owing [17] (solid curve) and additionally with energy- mom- 
entum conservation (dash-dotted) curve to data on J/ 1\) in 
d+Au collisions at yfs = 200 GeV [38 . The latter effect 
becomes important already at y — 2. Note, that the same 
mechanism will also affect the forward yields of charged 
hadrons at even smaller rapidities. 

This constitutes the baseline in the search for the ori- 
gin of anomalous suppression in high-energy AA collisions. 
Also, we also plot the predicted J/ip suppression in p+Pb 
collisions at LHC (see figure [TJ. 



3 Comover absorption and recombination of 
charmonium 

In the following section we will briefly describe the so- 
called co-mover model for the case of charmonium sec- 
ondary interactions. It was originally formulated to ex- 
plain the anomalous J /ip suppression in AA collisions at 
the SPS ^8 l [L9 l [20 l [2T]. Recently, predictions were made 
for RHIC including only comover dissociation [22] . 

As previously discussed, the possibility of recombina- 
tion of cc pairs in the dense medium cannot be neglected 
at RHIC, due to the considerable density of open charm 
at mid-rapidity. This recombination mechanism was first 
introduced in |23j neglecting spatial dependencies in the 
J/tjj and open charm densities. The approach has recently 
been improved in [241125) and charmonium recombination 
has also been considered in statistical hadronization mod- 
els with charm conservation [26,27]. The dynamics of char- 
monium recombination and dissociation has also been im- 
plemented in various scenarios in Monte-Carlo models (see 
e.g. HE]). 

The recombination mechanism has recently also been 
included in the co-movers interaction model (CIM) [29 . 
Assuming a pure longitudinal expansion and boost invari- 
ance of the system, the rate equation which includes both 
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Fig. 2. Results for J ftp suppression in Au+Au collisions at 
y/s = 200 GeV at mid- (upper figure) and forward (lower fig- 
ure) rapidities. Data are taken from [I]. 



dissociation and recombination effects for the density of 
charmonium at a given production point at impact pa- 
rameter s reads 



_dN J/4 ,(b,s,y) 



dr 



N c °(b, S ,y)N J/ ^b, S ,y) (9) 
- N c (b,s,y)No(b,s,y)\ , 



where N co , Nj/^ and N c ^ is the density of comovers, J/tj) 
and open charm, respectively, and <j C o is the interaction 
cross section for both dissociation of charmonium with 
co-movers and regeneration of J/ip from cc pairs in the 
system averaged over the momentum distribution of the 
participants. It is the constant of proportionality for both 
the dissociation and recombination terms due to detailed 
balance. 

Note, that the quantities that enter the rate equation, 
eq. ([9]), are densities of co- moving matter at the same im- 
pact parameter as the J/ip produced in the initial hard 
scattering. This is quite different from other approaches, 
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e.g. [23] and also the statistical hadronization models 26, 
I27j . where cc pairs in an extensive volume are allowed to 
recombine. In our case, the system is driven to a local 
equilibrium, given by 



N c (b,s,y)N s (b,s,y) 
N c °{b,s,y) 



(10) 



Equation © cannot be solved analytically. The sup- 
pression factor, i.e. the density of ratio of J/ip after the 
full evolution of the medium to the one at some formation 
time To, can be approximated by 

S co (b,s,y) = exp {- a co [N co (b,s,y) 

-C(y)N bln (b, S )S shad (b,s,y)] (11) 



x In [ 



N c 



N pp (0) 



}• 



where 



C(y) 



aXPda^/dy 



(12) 



PV 



Details of the model can be found in [25J . The quantities in 
eq. (|12[) are all related to pp collisions at the correspond- 
ing energy and are taken from experiment. Formulated as 
above, the extension of CIM with inclusion of recombina- 
tion effects does not involve any additional parameters. 

With a co = 0.65 nib 19,20,21 fixed from experiments 
at lower energies, we have calculated the suppression of 
J/ip in Cu+Cu and Au+Au collisions at ^/s = 200 GeV 
in [25]. We present results for the Au+Au case for both 
mid- and forward rapidities in fig. [2] Both the centrality 
and rapidity dependence are in good agreement with the 
data. 

The latter is due to a combination of stronger initial 
state effects at forward rapidities and the fact that re- 
combination effects are weaker than at mid-rapidity. In a 
natural way, the CIM explains both of the observed J/tp 
puzzles in A A collisions discussed earlier. 



4 J/i[) suppression at LHC 

Predictions for LHC energy of ^/s = 5.5 TeV can be 
readily done assuming that da C c/ dy\ y _ ~ 1 mb and 
<Jp p D = 59 mb which corresponds to C — 2.5; a co is kept 
constant. They are presented in fig. [3] 

These predictions are in stark contrast with predic- 
tions of the statistical hadronization model [30j . which 
favours a strong enhancement of J/ip for the most central 
collisions. 

The reason of the discreapancy is the local form of our 
rate equation: only comovers and open charm produced at 
the same impact parameter as the initial J/ip are allowed 
to interact. On the other hand, models assuming global 
equilibrium of the produced charm with the medium al- 
low for recombination of cc pairs from the whole volume 
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Fig. 3. Prediction for J/ip suppression in Pb+Pb collisions 
at mid-rapidity at LHC for different values of the parameter 
C. The upper line is the suppression due to gluon shadowing. 
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Fig. 4. Suppression of X at RHIC and LHC. 

of the fireball. A recent analysis suggests that the ther- 
mal relaxation times of charmed quarks are rather long at 
RHIC, t c ~ 5 — 7 fm/c [3T], but the situation is poorly 
known at LHC. Another point of dispute is the impor- 
tance of initial state effects, which amount to almost 50% 
of the suppression in the CIM. 

The large lever arm in energy and denisities of pro- 
duced will hopefully help us to understand and disentagle 
these issues in the future. 



5 Suppression of bottonium at LHC 

Finally, we present some results for bottonium at RHIC 
and at LHC, where excellent detector capabilities will al- 
low us to look in detail at the T family. Since the T 
state is much smaller than charmonium, it may provide a 
clearer probe of a thermalized system due to its presum- 
ably weaker interaction with matter. On the other hand, 
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Fig. 5. Centrality dependence of T suppression due to initial 
state effects (gluon shadowing) in Pb+Pb collisions at ^/s = 
5.5 TeV. 



due to strong binding energy, it may not dissolve at tem- 
peratures typical for heavy- ion collisions [32) . Predictions 
for LHC has previously been presented in [33J . 

Here, we would only like to discuss the impact of initial 
state effects on bottonium production (for a recent dis- 
cussion of final state medium effects we refer to [34]). The 
absorptive cross section for T is 40-50% smaller than the 
corresponding cross section for J ftp and ip' . Furthermore, 
energy-momentum conservation mechanisms are pushed 
to higher xp due to the large mass of the bottonium (cor- 
responding to 7 ~ 3 in eq. (O). We expect therefore neg- 
ligible nuclear absorption for RHfC and LHC kinematics. 

Shadowing for bottonium in pA collisions is shown in 
fig. [4] as a function of rapidity at RHIC and LHC energies. 
Since our model [17) does not contain antishadowing, the 
predictions for backward and mid-rapidities are uncertain 
up to the 10% level, but we expect shadowing at forward 
rapidity at RHIC. At LHC, a large suppression is pre- 
dicted for p+Pb collisions in most of the kinematics. 

The suppression of bottonium due to gluon shadowing 
in Pb+Pb collisions at LHC is shown in fig. [5] for sev- 
eral rapidities. The suppression is about 50% from mid- 
central to central collisions, and would be the same for 
all members of the T family. This establishes the base- 
line for further calculations of bottonium dissociation and 
recombination in the final state. 
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